{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "c95e1c49",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "等效解（辅助面）:\n",
      "87.40, 85.17, 75.78\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAScAAAEYCAYAAAAedjA5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA5WUlEQVR4nO3deVhUZf8/8PeZGWaGfRFBUBDFLfcFIdxIQ+nRNNukKKUey0pa1HKvqEwh08IU89IyS33ELXfDBcUdFwRFQVxQcQNFYUD2mbl/f/hlfpKAzMyZObN8Xtd1X+rhLG9gzsczZ+5z3xxjjIEQQkyMSOgAhBBSFypOhBCTRMWJEGKSqDgRQkwSFSdCiEmi4kQIMUlUnAghJomKEyHEJFFxIoSYJCpOhBCTpHVxOnjwIIYPHw5vb29wHIfNmzc/dZvk5GT07NkTMpkMbdq0wYoVK3SISgixJloXp9LSUnTr1g3x8fGNWv/q1asYNmwYBg4ciPT0dEyYMAHvvfcedu3apXVYQoj14PR58JfjOGzatAkjR46sd52pU6dix44dOHfunGbZG2+8gaKiIiQmJup6aEKIhZMY+gDHjh1DaGhorWVhYWGYMGFCvdtUVlaisrJS82+1Wo0HDx6gSZMm4DjOUFEJITpgjKGkpATe3t4Qifi7jW3w4pSXlwdPT89ayzw9PVFcXIzy8nLY2to+sU1MTAy+/fZbQ0cjhPDoxo0baNGiBW/7M3hx0sX06dMxadIkzb8VCgV8fX1x48YNODk5CZiMEPJvxcXF8PHxgaOjI6/7NXhxatasGfLz82sty8/Ph5OTU51XTQAgk8kgk8meWO7k5ETFiRATxfctF4P3cwoODkZSUlKtZXv27EFwcLChD00IMWNaF6eHDx8iPT0d6enpAB51FUhPT0dubi6AR2/JxowZo1n/ww8/RE5ODqZMmYILFy5g8eLFWLduHSZOnMjPd0AIsUxMS/v372cAnmiRkZGMMcYiIyNZSEjIE9t0796dSaVS1rp1a/bHH39odUyFQsEAMIVCoW1cQoiBGer81Kufk7EUFxfD2dkZCoWC7jkRYmIMdX7Ss3WEEJNExYkQYpKoOBFCTJJJdsIk5qO4uBjXrl3DjRs3cPv2bdy5cwf37t1DQUEBCgsLoVAoUFpairKyMlRUVECpVEKlUoHjOIhEItjY2EAmk8HOzg4ODg5wdnaGq6sr3N3d0bRpU3h5eaF58+bw8fGBn58fHBwchP6WiZFQcSJPxRjDrVu3cObMGZw7dw6ZmZm4ePEiLl26hPv372vW4zgOTZs2haenJ5o0aQI3Nzd4e3vDwcEBdnZ2kMlkkEgkEIvFAACVSgWlUomKigqUlZXh4cOHUCgUuH37NjIyMpCfn4+CggI8/plN06ZN0bZtW7Rv3x4dO3ZE586d0a1bNzRr1oyeu7QwVJzIExQKBVJSUnDkyBGsW7cO2dnZT93GwcEBHh4esLGxgUgkQlVVFcrKyiCXy+Hg4ABHR0d4eXnB19cXfn5+aNmyJWxsbJ663+rqauTl5SE3NxfXrl3DlStXcOnSJWRkZGDt2rUoKysD8Oh5zd69eyMwMBDBwcF49tln6SrLzFFXAoKysjJs2LABs2fPxsWLFxtcVyqVQqVSQaVSNWrfYrEYYrEYarUaSqVSs1wkEsHHxwddu3ZFt27d0L17dwQEBMDX17fRV0BqtRpXr17F2bNncfr0aZw8eRLHjx9HUVERxGIxevbsiUGDBuH5559H//79IZfLG7Vfoh1DnZ9UnKzU8ePH8eGHH2p6+guF4zhIJBJUV1cDAJo0aYJ+/fphwIABGDhwILp166bVMBxqtRoXLlzA4cOHkZycjP379yMvLw+2trYYOHAgXnzxRYwYMQLNmzc31Ldkdag4UXHSW0JCAsaNG4eSkhKhozRILBaDMQa1Wg1nZ2e88MILGDZsGIYOHYomTZpotS/GGM6fP4/ExETs3LkThw4dglKpRGBgIF577TW8/vrr8PPzM8w3YiUMdn7y2t/cQOjxFd0dO3aMeXh41PnIkbk0iUTCADCRSMQGDBjAFi9ezPLz83X6eRQWFrJVq1axl19+mcnlcgaA9e3bly1ZsoQVFhby+8O3EoY6P6k4WaCqqioWHh4ueFExRBOJRIzjOCYSidiQIUPY//73P1ZWVqbTz6m4uJitWrWK/ec//2EikYjJZDIWERHBkpOTmVqt5vm3YrmoOFFxeqqMjAzm4uIieAExVhOLxQwAc3R0ZJ9++inLzMzU+Wd3+/Zt9sMPP7C2bdsyAKxjx45s0aJFrKSkhMffkGWi4kTFqV5r164VvFAI3Wre+j333HNs69atTKVS6fSzVKvVLCkpib366qtMJBIxZ2dnNmXKFHbz5k2ef2uWg4oTFacnLFiwQPCiYGqt5mrK39+fLVu2jFVUVOj8871+/TqbPHkyc3JyYjY2Nuy9995jly5d4vE3aBmoOFFx0vjll18ELwKm3jiOYwBYs2bNWHx8vF5FSqFQsB9++IF5enoysVjMxowZwy5fvszjb9S8UXGi4sQ2bNgg+Elvjo3jOObl5cV+++03Vl1drfPPv6ysjP3yyy/My8uLSSQS9sEHH7Dbt2/z+Bs2T1ScrLg4nT17VvAT3NxbzZVU27Zt2bZt2/T6NK6srIzNmzePubm5MTs7O/b1119b9Y1zKk5WWJzKy8uZn5+f4Ce2JbWae1KDBg1i586d0+v3U1hYyKZOncpkMhnz8vJiK1eutMouCFScrKw4TZ8+XfAT2ZKbRCJhIpGITZo0iRUXF+v1u7p69Sp77bXXGADWr18/dvbsWZ5eBeaBipOVFKeLFy8KfuJaUxOJRKxZs2Zs8+bNev/u9u7dyzp06MAkEgmbNm2azp1DzQ0VJysoTjX/+1IzfoECwF5//XV29+5dvX6HlZWVbNasWUwqlbI2bdqwgwcP8vTqMF1UnCy4OF2+fFnwE5Tao/tRrq6ubOPGjXr/TrOyslifPn0Yx3FswoQJFn0VRcXJQovThAkTBD8pqf3/VvOpXmRkpN73opRKJZs/fz6TyWTsmWeeYWlpafy8aEwMFScLK04VFRVMKpUKfjJSq7uJxWLWsmVLdurUKb1/1+fPn2fdunVjUqmUxcXFWdwneoY6P2n2FQEcPHgQcrkcVVVVQkch9VCpVLh58yaeffZZLFy4sNY45trq2LEjjh8/jvHjx2PChAkYOXIkCgsLeUxroXgtdQZiSVdOn376qeBXBdS0b+Hh4bx0tNy6dStzcXFhrVq1YqdPn+bhFSU8unKyAH5+fvjll1+EjkF0sH79egQGBuLKlSt67Wf48OE4ffo0XF1d0adPH6xevZqnhJaHipMRPHz4EBzH4fr160JHITpSq9W4dOkSevXqhf379+u1r1atWuHw4cMYNWoU3n77bUybNg1qtZqnpJaDipOBnT9/Ho6OjkLHIDxQKpUoKSnB4MGD8fvvv+u1L1tbW6xYsQLz58/H3Llz8eqrr6K0tJSnpBaC1zeJBmKu95y2bNki+L0SaoZrM2fO5OWTt61btzJ7e3vWu3dvncdGFxLdczIzcXFxeOmll4SOQQxo9uzZeOeddzTTWulq+PDhOHjwIHJzc9GnTx/k5OTwlNC8UXEygMmTJ2PixIlCxyBGsHLlSowcORLl5eV67adnz544duwYOI5D3759kZGRwVNC80XFiWdjx47FvHnzhI5BjIQxhsTERLzwwgt6zwfYqlUrHDlyBF5eXggJCcHJkyd5SmmeqDjxKCIiAsuXLxc6BjEytVqNI0eO4Pnnn4dCodBrXx4eHti3bx86dOiA0NBQHDt2jKeU5oeKE0/efvttrFmzRugYRCAqlQqnT5/GoEGDUFRUpNe+XFxcsGvXLnTt2hVDhgxBSkoKPyHNjE7FKT4+Hn5+fpDL5QgKCsKJEycaXD8uLg7t27eHra0tfHx8MHHiRFRUVOgU2BSNGzeOOtMRqFQqnDlzBoMHD0ZxcbFe+3J0dMQ///yD7t27IywsDKmpqTylNCPafryXkJDApFIpW758OTt//jx7//33mYuLS70fga5evZrJZDK2evVqdvXqVbZr1y7m5eXFJk6c2OhjmnJXAhqxktq/m1gsZn379mWlpaV6v76Ki4tZUFAQa9Kkid7DChuKyYxKEBgYyKKiojT/VqlUzNvbm8XExNS5flRUFBs0aFCtZZMmTWJ9+/Zt9DFNtTgtXrxY8BOBmmk2sVjMhgwZwiorK/V+nT148IB17dqVNW/enF27do2HVy6/TKKfU1VVFVJTUxEaGqpZJhKJGrxx16dPH6Smpmre+uXk5GDnzp0YOnRovceprKxEcXFxrWZqdu3ahfHjxwsdg5golUqFvXv3IjIyUu9HU1xdXbFr1y7IZDKEhYXh/v37PKU0cdpUslu3bjEA7OjRo7WWT548mQUGBta73YIFC5iNjY1myugPP/ywweNER0fX+b+RqVw50ciV1LRpn3/+OS+vu0uXLjF3d3fWt29fVl5ezss++WASV066SE5Oxpw5c7B48WKcPn0af//9N3bs2IFZs2bVu8306dOhUCg07caNG4aO2Wjl5eVo06aN0DGIGZk/fz4WLVqk937atGmD7du3IzU1Ff/973/1GmPKLGhTySorK5lYLGabNm2qtXzMmDFsxIgRdW7Tr18/9sUXX9RatnLlSmZra8tUKlWjjmtK95xgAv8TUzO/xnEc2759Oy+vwXXr1jEAbNasWbzsT18mceUklUrRq1cvJCUlaZap1WokJSUhODi4zm3KysogEtU+jFgsBgCzq/xdu3YVOgIxY6+//jrOnTvHy36+/fZbfPXVV9iyZQsPyUyUttUsISGByWQytmLFCpaZmcnGjRvHXFxcWF5eHmOMsdGjR7Np06Zp1o+OjmaOjo5szZo1LCcnh+3evZv5+/uzUaNGNfqYpnDlVN99MGrUGtskEgnz8fFhBQUFer8e1Wo1e+WVV5ijoyPLysri4RWuO5PpSsAYYwsXLmS+vr5MKpWywMBAlpKSovlaSEgIi4yM1Py7urqaffPNN8zf35/J5XLm4+PDxo8fzwoLCxt9PKGLU1pamuAvbGqW0cRiMRs0aBBTKpV6vy6Li4tZhw4dWKdOndjDhw95eKXrxqSKk7EJWZyUSqXgL2hqltU4jmPTp0/n5fV5/vx5Zm9vz8aMGcPL/nRhEvecrJGXl5fQEYiFYYwhJiYG27dv13tfHTt2xJIlS/DXX3/hr7/+4iGd6aDi1ICYmBjcu3dP6BjEAnEch7fffpuXbjJvv/02IiMjERUVhcuXL/OQzjRwjJn+R2bFxcVwdnaGQqGAk5OTUY559+5deHp6GuVYxDpJJBL07t0bBw8ehEQi0WtfJSUl6NGjB5o2bYpDhw7pvT9tGOr8pCunelBhIoamVCqRkpKCOXPm6L0vR0dH/PXXXzhx4gTmzp3LQzoTwOsdLAMx9g3xiRMnCn7TlJr1NJFIVOsTb31MmzaN2djYsIyMDF721xiGOj/pbd2/FBQUoGnTpgY9BiGPE4vFaNmyJTIyMmBnZ6fXviorK9GzZ0/Y29vj2LFjmg7PhkRv64zEw8ND6AjEyqhUKly7dg0zZszQe18ymQy//fYbTp06xcvzfEKi4vSYZcuWmd0jNcQyqNVq/PLLLzhy5Ije+woODsb48ePx5Zdf4ubNmzykEwa9rfs/arXaKJfAhNRHLBajVatWOHfuHGQymV77KioqQocOHTBgwACsW7eOp4R1o7d1Bvb4AHqECEGlUuHKlSuIiYnRe18uLi748ccfsX79euzbt4+HdMZHV04A7t27R/eaiMmQSCQ4f/482rVrp9d+GGPo168fSkpKkJaWZrB3BnTlZED+/v5CRyCklvHjx+t9/5PjOPz888/IyMgwy/kUrb44nT59Wu+ZWgnhk1KpRFJSEjZt2qT3vgIDAxEREYHo6GiUlpbykM54rL449erVS+gIhDxBJBLhs88+42V+x++//x4FBQVYsGABD8mMx6qL0969e4WOQEid1Go1bt26hbi4OL331apVK3zwwQf48ccf9Z6N2JisujgNHjxY6AiE1IsxhlmzZuHu3bt672vmzJmoqKjgpdgZi9UWp/379wsdgZCnqqysbHCmosZq1qwZPvroI8TFxUGhUPCQzPCstjgNGjRI6AiEPJVKpcKSJUuQk5Oj974mT56M8vJyLF68mIdkhmeVxYmPGTAIMaZvvvlG7314eXnhv//9L37++WeUl5frH8rArLI41TeNFSGmSKlUYtWqVcjOztZ7X59//jkKCgqwcuVKHpIZltUVp4cPH+Lhw4dCxyBEK2KxGN99953e+2nTpg1GjhyJuLg4k3/I3eqKU1hYmNARCNGaUqnEmjVreBkjfMKECcjKysKePXt4SGY4Vlecjh49KnQEQnQiFosRGxur93769++Prl27mvyNcasqTuY++BaxbkqlEn/++Sdu376t1344jsNHH32Ebdu28TL7i6FYVXH65JNPhI5AiF4YY7z8JxsREQG5XI4//viDh1SGYTXFyVw6nhHSEJVKhfj4eL0f4nVyckJ4eDhWrFgBtVrNUzp+WU1xGjVqlNARCOFFSUkJL7P7vvPOO7h69SoOHz7MQyr+WU1x2r17t9ARCOHNggUL9O4K0L9/f/j5+WHVqlU8peKXVRSnS5cuCR2BEN4wxpCdnY0DBw7otR+O4xAREYH169ejqqqKp3T8sYri9OabbwodgRBeSSQSLFmyRO/9vPHGGygqKjLJ4YOsojilpqYKHYEQXimVSmzcuBEFBQV67adz585o3749NmzYwFMy/lh8cdK3TwghpkqlUul9v4jjOLz88svYunUrlEolT8n4YfHF6f333xc6AiEGwRjD77//rvd+Ro4cifv37yMlJYWHVPyx+OK0c+dOoSMQYjDnzp1DRkaGXvvo3bs3mjZtiu3bt/OUih8WXZxMtXMZIXyRSCR6v7UTiUR44YUXkJiYyFMqfuhUnOLj4+Hn5we5XI6goCCcOHGiwfWLiooQFRUFLy8vyGQytGvXzihXNKbcNZ8QPtSM9aRvn6ewsDCcOXMG+fn5PCXjAdNSQkICk0qlbPny5ez8+fPs/fffZy4uLiw/P7/O9SsrK1lAQAAbOnQoO3z4MLt69SpLTk5m6enpjT6mQqFgAJhCodAqa4sWLRgAatQsvh09elSrc+Pf7ty5wwCwNWvWaL2trufn02hdnAIDA1lUVJTm3yqVinl7e7OYmJg61//1119Z69atWVVVlc4hdf3mhX7BUKNmjCaRSNikSZN0Pr9qPPPMM2zcuHFab2eo4qTV27qqqiqkpqYiNDRUs0wkEiE0NBTHjh2rc5utW7ciODgYUVFR8PT0ROfOnTFnzhyoVCptDq01ut9ErIVSqcT69ev1fms3YMAAHDp0iKdU+tOqOBUUFEClUsHT07PWck9PT+Tl5dW5TU5ODjZs2ACVSoWdO3fiq6++wvz58/H999/Xe5zKykoUFxfXatpau3at1tsQYq5u3Lih98Qd/fr1Q1ZWFu7fv89TKv0Y/NM6tVoNDw8PLF26FL169UJ4eDhmzpzZYNf7mJgYODs7a5qPj4/Wx503b54+sQkxK2KxGNu2bdNrHzUTf5hKfyetipO7uzvEYvETd/Tz8/PRrFmzOrfx8vJCu3btIBaLNcueeeYZ5OXl1fuw4fTp06FQKDRNl9H6Tp8+rfU2hJgrlUqFLVu26LWP1q1bw93dHSdPnuQplX60Kk5SqRS9evVCUlKSZplarUZSUlK90y317dsXly9frnUP6OLFi/Dy8oJUKq1zG5lMBicnp1qNENKwkydP4sGDBzpvz3EcAgICcOrUKR5T6UHbO+gJCQlMJpOxFStWsMzMTDZu3Djm4uLC8vLyGGOMjR49mk2bNk2zfm5uLnN0dGQff/wxy87OZtu3b2ceHh7s+++/b/Qxtf00oKqqSvBPUKhRE6KtX79euxP6X6ZPn86aN2+u1TYm05WAMcYWLlzIfH19mVQqZYGBgSwlJUXztZCQEBYZGVlr/aNHj7KgoCAmk8lY69at2ezZs5lSqWz08bT95letWiX4i4QaNWM3Gxsb9sEHHzT6vKpLQkICA8AKCgoavY1JFSdj0/abHzJkiOAvFGrUhGh+fn56nWvnzp1jANiBAwcavY1J9HMyF6Y6JjIhhnbt2jW9hglq27YtxGIxMjMzeUylG4ssTmVlZUJHIEQwBw8e1HlbqVSK1q1b4+LFizwm0o1FFidCrJWNjY3es1q3bdvWJMbdp+JEiAWprq7W+xGUVq1a4erVqzwl0p3FFSeaPJNYu3PnzqGiokLn7f38/HD9+nW9n9XTl8UVJ5qfjlg7pVKJs2fP6ry9j48PHj58KPh/9BZXnEzpqWpChMBxnF6PbzVv3hyA8JODWFxx0nc8ZULMnUQiQVpams7b1zwnW99II8ZiccXJFG7kESKk6upqvYqTh4cHAOg9J56+LK44CV3tCTEF58+f1/mGtqOjIyQSCRUnvlVWVgodgRDBlZWV4ebNmzpty3EcXFxcUFRUxG8oLVlccSKEPJKdna3TdsXFxSgoKBB8NFkqToRYII7jtHoERalUIjU1FT/99BPat28PAHp1R+CDRNCjE0IMwsbGBleuXKnza2q1GtevX0dmZibS09Nx4MABHDlyBGVlZRCJTOd6hYoTIRaouroa6enpSEtLQ05ODi5cuICsrCycPXsWFy9e1NybFYvFYIxpRqo1pVmLqDgRYmFsbGxQXV2Nffv2oWfPngAe9X0CHr19e5yhp2jTBxUnQsyIWCyGWCyGUqmsdZUjFovRvHlztGnTBv7+/li1ahXKy8s1X/93UTIHVJwIMQESiQQikQgqleqJqxmpVAoPDw/4+vqiZcuWaNGixRPN09Oz1gxHzz77LMaOHWvsb4NXVJwIMQCRSASJRALGGJRK5RMdIkUiEVxcXODp6QkvLy94eXnB09MTzZo1g5eXF7y9vTXLnZ2dwXGcVsevb6o2c0LFiZCnqLmqYYxBpVLVedNYLBbDxcUF7u7u8PT0hKenJ9zd3eHu7g4PDw94eHigadOmmj+bNGlS60qHb/+eldscUXEiFo/jOEgkEnAc12CBAR4VGWdnZ7i4uKBJkyaaAuPq6go3Nze4ubmhSZMmaNKkCdzc3ODu7o4mTZrA0dFR66sbQ3J3dxc6gt6oOBGTJRKJIBaLNSe9Wq2GSqVq8JkxGxsbODo6wsnJCS4uLnBzc4OLi4um4Li6umr+XvPvx/+0t7c3qSKjKzc3N6Ej6M3iitPgwYOxZ88e2Nraory8XPPpRn1TnxP91VyV1JzUNf1mGvsxtVwuh52dHezt7TWFxdnZGU5OTnBwcICTkxMcHR017d+zQT/eZDKZIb9Vs+Hg4ACRSGRS/Za0ZXHFKSgoCFlZWbhw4QJ2796Ny5cv4+rVq8jJycGlS5dw48YNVFdX17mtSCSCSCR64kSreSsg9LClNWry1fTmfTxvjZrcNYVC1+w2NjaQyWSQy+WQyWSws7PTFJKaZmdnB1tbW83f6/q6vb09HBwc6vzTlHolWwqO42Bvb4+SkhKd9+Hk5MRjIu1ZXHFyc3NDYWEh7O3t8fLLLz/xdcYYSktLUVRUhMLCQhQVFdX6e1lZGcrLy1FRUYHy8vJaf6+oqIBSqUR1dXWtplQqUVVVVasQ/PvPxwvI438Xi8Wagljzd7FYDIlEovmzptnY2MDGxqbev9c0mUym+VMqlUIqlWr+XlNkappcLtcsq/m7XC6Hra0t5HI5FQ4zZmdnp1dxEvqmusUVJ3d3d5SWlqK8vBy2trZPfJ3jODg4OMDBwQEtWrQQICEhxmFnZ6fX9n5+fvwE0ZHF/bdYM4rf3bt3BU5CiLDkcrle23fu3JmnJLqxuOLk5eUFALhz547ASQgRlo2NjV7bBwcH85RENxZXnGpmjtB1FEBCLIW+nTwHDx7MUxLdWFxxcnNzg729PXJzc4WOQohZc3FxEfT4FlecOI6Dn58fzcJCrJ6pdH3RlcUVJwDw9/fH5cuXhY5BiKDMcZiUx1lkcWrTpo1W4ycTYon0eSrCkA8lN5ZFFqdnnnkG165dQ0VFhdBRCBHM44PNaSsgIIDHJLqxyOLUsWNHqNVqXLhwQegohAhGn+L0+uuv85hENxZZnDp16gQAyMjIEDgJIcJ5+PChztuawiiaOhWn+Ph4+Pn5QS6XIygoCCdOnGjUdgkJCeA4DiNHjtTlsI3m7OyMVq1a4cyZMwY9DiGmSqlU6nVbQ+huBIAOxWnt2rWYNGkSoqOjcfr0aXTr1g1hYWFPfVzk2rVr+OKLL9C/f3+dw2qjZ8+eSE1NNcqxCDE1hYWFQkfQm9bF6aeffsL777+Pd999Fx07dsSSJUtgZ2eH5cuX17uNSqXCW2+9hW+//RatW7fWK3BjBQQE4NSpUyY99Q0hhlJQUKDztjVPWQhNq+JUVVWF1NRUhIaG/v8diEQIDQ3FsWPH6t3uu+++g4eHR6Pfx1ZWVqK4uLhW01ZQUBAePnyIzMxMrbclxNzp8+C7KdxvArQsTgUFBVCpVE+M8+Lp6Ym8vLw6tzl8+DB+//13LFu2rNHHiYmJgbOzs6b5+PhoExMAEBgYCLFYjKNHj2q9LSHmrr7zsTEmT57MYxLdGfTTupKSEowePRrLli3TasD16dOnQ6FQaNqNGze0Pra9vT169uyJQ4cOab0tIebu1q1bOg8U6ODgwHMa3Wg12Jy7uzvEYjHy8/NrLc/Pz69znqwrV67g2rVrGD58uGZZzZjGEokE2dnZ8Pf3f2K7mlEa9RUSEoKEhATNSJSEWIubN29CIpFo3UvcFHqG19CqtEqlUvTq1QtJSUmaZWq1GklJSXWO/dKhQwdkZGQgPT1d00aMGIGBAwciPT1dp7dr2hg4cCBu3ryJS5cuGfQ4hJiaa9eu6fRs3bhx4wyQRjdaD9M7adIkREZGIiAgAIGBgYiLi0NpaSneffddAMCYMWPQvHlzxMTEQC6XPzGaXk3/CWOMsjdgwADY2Nhgz549aNeuncGPR4ipuHjxok4zr8TExBggjW60Lk7h4eG4d+8evv76a+Tl5aF79+5ITEzU3CTPzc01mUHxHRwc0LdvXyQmJiIqKkroOIQYBWMMOTk5Om3r7OzMcxrdccwMBn0pLi6Gs7MzFAqF1tPV/Pjjj4iOjsb9+/frnPCAEEtz69YtnSbv8PX1xfXr17XeTp/zsyGmcYljQC+++CLKy8tr3ScjxJJlZWXptN2iRYt4TqIfiy9OHTp0QNu2bbF582ahoxBiFOfOndPp1srjn6qbAosvThzH4dVXX8XmzZvNfmRAQhrj7NmzWncJMMVp3C2+OAHAa6+9hvv372Pfvn1CRyHE4E6ePInq6mqttjG1t3SAlRSnnj17ok2bNkhISBA6CiEGVVFRodM9p/fee88AafRjFcWJ4zi8+eab2Lhxo16jAxJi6tLS0rQeicOUug88ziqKEwCMHj0axcXF2LJli9BRCDGYlJQUrW+Gr1u3zkBp9GM1xalt27bo27dvg+NOEWLuDh06pPVzpEOGDDFQGv1YTXECHo1Ts3fvXppwk1gkxhiSk5O1elsXEhJiwET6sariNGrUKDg5OWHp0qVCRyGEd5mZmVoPz7t161YDpdGfVRUne3t7REZG4rfffqM57YjFSUpK0vp+E5+Pm/DNqooTAERFRaGgoIC6FRCLk5iYqNX9pu3btxswjf4s/sHfugwbNgw3btzAmTNnaBA6YhHKy8vh6uqKysrKRm/D16lPD/7y6PPPP0dGRgb27NkjdBRCeLFv3z6tCtOMGTMMmIYfVlmcBg4ciF69eiE2NlboKITwYuvWrbCxsWn0+rNnzzZgGn5YZXHiOA4zZszA/v37G5zSihBzoFKpsGHDhkY/T/fGG28YOBE/rPKeE/Bo7POuXbuiRYsWSExM5GWfhAghOTkZAwcObPT6KpWK19Fq6Z4Tz0QiEb766ivs2rWL5rYjZi0hIQESSeNG3H755ZdNZhjtp7HaKyfg0dVTjx494Orqiv3799Mnd8TsVFZWwsPDo9GzYvN91QTQlZNBiEQizJ49GwcOHMCuXbuEjkOI1nbs2NHowvTee++ZzVUTYOVXTsCjvh4hISEoKipCWlqaSU0qSMjTDBs2DLt3727UKK+GOtXpyslAOI7D/PnzkZGRgd9//13oOIQ02s2bN/HPP/80qjDFxcUZPhDPrL44AUDv3r0xevRozJw5E0VFRULHIaRRli5d2ui3aZ999pmB0/CPitP/iY2NRUVFBaKjo4WOQshTVVVV4ddff23U8CgnT540QiL+UXH6P97e3oiOjsaiRYuQlpYmdBxCGrRu3ToUFBQ8db2WLVsiICDACIn4Z/U3xB9XXV2Nnj17Qi6XIyUlhW6OE5PEGEO3bt2QmZn51CunyspKSKVSg+ahG+JGYGNjg6VLlyI1NRW//PKL0HEIqdOePXuQkZHx1MIUExNj8MJkSHTlVIdPP/0Uv/32G86ePYs2bdoY/HiEaGPAgAE4duzYUz+lM9apTVdORjRnzhx4eXnh3Xff1XqaHUIM6cCBAzh06NBTC9O9e/eMlMhwqDjVwcHBAX/88QeOHDmC+fPnCx2HEACProS+/PLLpz5H9+WXX8Ld3d1IqQyH3tY1YMqUKYiLi8Px48fRo0cPox2XkLokJibiP//5T4PrSCQSraci15ehzk8qTg2orKxEcHAwSktLkZqaCgcHB6Mdm5DHqVQqdO/eHVlZWQ3eaigrK4Otra0Rk9E9J0HIZDIkJCTg1q1bGD9+vNFuMBLyb3/99RfOnTvXYGHauHGj0QuTIVFxeop27dph6dKlWLlyJZYtWyZ0HGKFSkpKMHXq1AaH9BkyZAheeeUVI6YyPCpOjRAREYEPP/wQn3zyidk+CkDM13fffYf79+83eOVuiUP+6FSc4uPj4efnB7lcjqCgIJw4caLedZctW4b+/fvD1dUVrq6uCA0NbXB9UxUXF4cePXrglVdeQV5entBxiJXIyMjAzz//DLVaXe86Dx8+NGIi49G6OK1duxaTJk1CdHQ0Tp8+jW7duiEsLAx3796tc/3k5GS8+eabmskEfHx8MGTIENy6dUvv8MYkk8mwceNGqFQqvPLKKzRjMDE4lUqFsWPHNvh27tixY7C3tzdiKiNiWgoMDGRRUVGaf6tUKubt7c1iYmIatb1SqWSOjo7szz//bPQxFQoFA8AUCoW2cXmXkpLCZDIZe+utt5harRY6DrFgcXFxDEC9be7cuUJHZIwZ7vzU6sqpqqoKqampCA0N1SwTiUQIDQ1t9BRLZWVlqK6uhpubW73rVFZWori4uFYzFUFBQfjzzz+xevVqGl6FGMylS5cwderUer/+4osvYvLkyUZMZHxaFaeCggKoVCp4enrWWu7p6dno+zBTp06Ft7d3rQL3bzExMXB2dtY0Hx8fbWIaXHh4OGJjYzFr1iwsXbpU6DjEwiiVSkRERNTbbaBFixbYtm2bkVMZn1E/rYuNjUVCQgI2bdoEuVxe73rTp0+HQqHQtBs3bhgxZeNMmTIFUVFR+Oijj/D3338LHYdYkFmzZiE1NbXe5+dM8XwwhMZNdvV/3N3dIRaLkZ+fX2t5fn4+mjVr1uC28+bNQ2xsLPbu3YuuXbs2uK5MJoNMJtMmmtFxHIcFCxbg3r17ePPNN7F9+3YMHjxY6FjEzCUnJ2PWrFn1dhtozHjhlkKrKyepVIpevXohKSlJs0ytViMpKQnBwcH1bjd37lzMmjULiYmJZjsqX13EYjFWrlyJ0NBQvPTSSzh48KDQkYgZu3PnDkaNGlXvuOBlZWXWNQCitnfQExISmEwmYytWrGCZmZls3LhxzMXFheXl5THGGBs9ejSbNm2aZv3Y2FgmlUrZhg0b2J07dzStpKSk0cc0pU/r6lJWVsZCQ0OZvb09O3TokNBxiBmqrKxkwcHBTCKR1PnJXGFhodAR62Wo81Pr4sQYYwsXLmS+vr5MKpWywMBAlpKSovlaSEgIi4yM1Py7ZcuWdf6wo6OjG308Uy9OjDFWWlrKBg4cyOzt7dn+/fuFjkPMzEcffcQ4jqvzXLl3757Q8RpkqPOTRiXgUVlZGV566SUcPnwYmzZtwgsvvCB0JGIGFi5ciE8//bTOr92/f7/BbjemgEYlMAN2dnbYtm0bBg8ejBEjRmDdunVCRyImbtu2bfXOKadQKEy+MBkSFSeeyeVybNy4EeHh4XjjjTewaNEioSMRE5WSkoLXX3+9zq+Vl5eb9LsEY9CqKwFpHBsbG/z555/w9PTEJ598guvXr+OHH35o9OysxPKdO3cOYWFhUCqVT3QbUKlU9FoBFSeDEYlEmDdvHnx8fDBx4kTk5OTgr7/+styHNEmjZWdnY+DAgSgrK6vVC9zf3x+XL18WMJlpofJsYJ999hk2b96MXbt2oX///sjNzRU6EhFQdnY2BgwYgKKiolodKsPDw6kw/QsVJyMYMWIEjhw5ggcPHiAgIAAHDhwQOhIRwPnz59GvXz88ePCgVmFavHgxEhISBExmmqg4GUm3bt1w8uRJdO7cGc8//zx+/PFHGpPcihw/fhx9+/ZFYWFhrcJ05swZfPTRRwImM11UnIyoadOm2L17N7744gtMmTIFL730Eh48eCB0LGJgO3bswHPPPQeFQqG5xyQWi1FZWfnU50ytGRUnI5NIJIiNjcW2bdtw5MgRdO3aFcnJyULHIgayePFiDB8+vNbIqREREVAqlZBKpQImM31UnATy4osv4syZM2jTpg0GDRqEqVOnorKyUuhYhCfV1dX4+OOPERUVVevt+/79+7F69WoBk5kPKk4CatGiBZKSkhAbG4uff/4ZAQEBSE1NFToW0dOdO3cwcOBAxMfHa5Z5eHhAqVTiueeeEy6YmaHiJDCxWIwpU6bg1KlTkEgkCAoKwrRp01BeXi50NKKDpKQkeHt748iRI5pl8fHxyM/Pt67hTnhAxclEdO3aFSdOnMC3336Ln3/+GV26dMHu3buFjkUaqaqqCjNmzKg1/LSzszPKysowfvx4AZOZLypOJsTGxgYzZ87E2bNn4evri7CwMLz22mvUcdPEZWRkQCaTISYmRrNs9erVKCoqsqjpwY2NipMJat++PZKSkrB69WocPXoU7du3R3R0tMVOnmiuKisrMXbs2FrdAfr37w+VSoWIiAgBk1kGKk4miuM4REREIDs7GxMmTEBsbCzatm2LpUuXWtU40qZqy5YtkMvlWL58uWZZbm4uDh48SA/t8oR+iibO0dERMTExyM7OxqBBg/DBBx+gY8eOWLNmTb1TBxHDOX/+PDiOw8iRIzXLtm/fDsaYyU1hZu6oOJkJPz8/rF69GmlpaWjXrh0iIiLQpUsXrF69mq6kjODq1avgOA6dO3fWLPvpp5/AGMOwYcMETGa5qDiZme7du2P79u1ISUlBq1at8Pbbb6N9+/b49ddfUVZWJnQ8i5Oeng6O49C6dWvNsi+++AKMMUycOFHAZJaPipOZCgoKwo4dO5CamoqAgAB8/PHH8PX1xcyZM3Hz5k2h45m9pUuXguM49OjRQ7Ps888/B2MMP/74o4DJrAiv0yUYiDnMviK0K1eusM8++4w5ODgwsVjMXnnlFbZr1y6mUqmEjmY2SktLWfPmzZ+Y/SQ+Pl7oaCbNUOcnXTlZiNatWyMuLg63b9/GggULcPHiRYSFhaF169aIjo6mgcwakJWVhSlTpsDe3h63bt3SLD979iwYY9SJUiA0NZSFYowhJSUFy5cvx9q1a1FSUoKgoCC88cYbeO2119CiRQuhIwrq2rVr2LBhA9asWYPTp0/DxcUFRUVF6N+/P83crCVDnZ9UnKxAWVkZtm7dijVr1uCff/5BdXU1goKCMHLkSAwfPhwdO3YEx3FCxzQoxhjOnj2LrVu3YsuWLUhNTYVMJsOLL76IiIgIDBs2DDKZTOiYZomKExUnXhQVFWH79u34+++/sXv3bpSWlmoelQkNDcXAgQPRtGlToWPy4s6dO9i/fz/27NmD3bt34/bt23B0dMQLL7yAV199FUOHDoWjo6PQMc0eFScqTryrqKhAcnIyEhMTsXv3bmRlZQEAOnbsiH79+qFPnz4ICgpCu3btTL7Xs0qlwoULF3D8+HEcPXoUhw8fRnZ2NgCgS5cuGDJkCIYOHYp+/frRIG88o+JExcngbt26heTkZBw4cABHjhxBZmYmgEe91Hv06IHu3bujc+fO6NSpEzp06CDIbLSMMRQUFODChQvIzMxERkYG0tPTkZ6ejtLSUnAch06dOqFfv34ICQnBc889h2bNmhk9pzWh4kTFyeiKiopw6tQpnDp1CmlpaThz5gwuXboEtVoNAHB1dYW/vz/8/Pzg6+uL5s2bw8vLC02bNoW7uztcXV3h5OQER0dHSCQNT5FYXV2Nhw8fQqFQoLCwEAUFBcjPz0deXh5u3bqF3NxcXLt2DVeuXIFCoQDwaCysdu3aoWvXrujVqxd69eqFgIAAeo0YGRUnKk4moaKiApcuXcKFCxdw+fJlXLlyBdevX0dubi5u3bqF0tLSOrezsbGBTCaDjY2N5i2iSqWCUqlERUVFvY/gODo6onnz5vDx8UHLli3h7++Ptm3bokOHDmjTpg3dxDYBhjo/acZfohW5XI4uXbqgS5cuT3yNMYaHDx/i3r17uH//PgoLC1FcXIySkhJUVFRoipBSqQTHcRCJRJqiZWtrCycnJzg5OcHNzQ1ubm7w9PSEnZ2dAN8lMQVUnAhvOI6Do6MjHB0daz2LRoguTPsjGEKI1aLiRAgxSVScCCEmiYoTIcQk6VSc4uPj4efnB7lcjqCgIJw4caLB9devX48OHTpoPunZuXOnTmEJIdZD6+K0du1aTJo0CdHR0Th9+jS6deuGsLAw3L17t871jx49ijfffBNjx45FWloaRo4ciZEjR+LcuXN6hyeEWC6tO2EGBQWhd+/eWLRoEQBArVbDx8cHn3zyCaZNm/bE+uHh4SgtLcX27ds1y5599ll0794dS5YsadQxqRMmIabLJDphVlVVITU1FdOnT9csE4lECA0NxbFjx+rc5tixY5g0aVKtZWFhYdi8eXO9x6msrERlZaXm3zWPKxQXF2sTlxBiBDXnJd8Pm2hVnAoKCqBSqeDp6VlruaenJy5cuFDnNnl5eXWun5eXV+9xYmJi8O233z6xnKbeIcR03b9/H87OzrztzyR7iE+fPr3W1VZRURFatmyJ3NxcXr95QysuLoaPjw9u3Lhhdm9HzTW7ueYGzDe7QqGAr68v76NUaFWc3N3dIRaLkZ+fX2t5fn5+vcNSNGvWTKv1AUAmk9X5QKezs7NZ/dJq1DwzZo7MNbu55gbMNzvfY35ptTepVIpevXohKSlJs0ytViMpKQnBwcF1bhMcHFxrfQDYs2dPvesTQgigw9u6SZMmITIyEgEBAQgMDERcXBxKS0vx7rvvAgDGjBmD5s2bIyYmBgDw2WefISQkBPPnz8ewYcOQkJCAU6dOYenSpfx+J4QQi6J1cQoPD8e9e/fw9ddfIy8vD927d0diYqLmpndubm6ty7s+ffrgf//7H7788kvMmDEDbdu2xebNm2tN6/w0MpkM0dHRZjd2j7nmBsw3u7nmBsw3u6Fym8Vgc4QQ60PP1hFCTBIVJ0KISaLiRAgxSVScCCEmyWSKk7kOw6JN7mXLlqF///5wdXWFq6srQkNDn/p9GpK2P/MaCQkJ4DgOI0eONGzAemibu6ioCFFRUfDy8oJMJkO7du3M4vUCAHFxcWjfvj1sbW3h4+ODiRMnoqKiwkhpHzl48CCGDx8Ob29vcBzX4HOxNZKTk9GzZ0/IZDK0adMGK1as0P7AzAQkJCQwqVTKli9fzs6fP8/ef/995uLiwvLz8+tc/8iRI0wsFrO5c+eyzMxM9uWXXzIbGxuWkZFh0rkjIiJYfHw8S0tLY1lZWeydd95hzs7O7ObNm0bNzZj22WtcvXqVNW/enPXv35+99NJLxgn7GG1zV1ZWsoCAADZ06FB2+PBhdvXqVZacnMzS09ONnFz77KtXr2YymYytXr2aXb16le3atYt5eXmxiRMnGjX3zp072cyZM9nff//NALBNmzY1uH5OTg6zs7NjkyZNYpmZmWzhwoVMLBazxMRErY5rEsUpMDCQRUVFaf6tUqmYt7c3i4mJqXP9UaNGsWHDhtVaFhQUxD744AOD5vw3bXP/m1KpZI6OjuzPP/80VMR66ZJdqVSyPn36sN9++41FRkYKUpy0zf3rr7+y1q1bs6qqKmNFrJe22aOiotigQYNqLZs0aRLr27evQXM2pDHFacqUKaxTp061loWHh7OwsDCtjiX427qaYVhCQ0M1yxozDMvj6wOPhmGpb31D0CX3v5WVlaG6utro03rrmv27776Dh4cHxo4da4yYT9Al99atWxEcHIyoqCh4enqic+fOmDNnDlQqlbFiA9Ate58+fZCamqp565eTk4OdO3di6NChRsmsK77OT8FHJTDWMCx80yX3v02dOhXe3t5P/CINTZfshw8fxu+//4709HQjJKybLrlzcnKwb98+vPXWW9i5cycuX76M8ePHo7q6GtHR0caIDUC37BERESgoKEC/fv3AGINSqcSHH36IGTNmGCOyzuo7P4uLi1FeXg5bW9tG7UfwKydrFRsbi4SEBGzatAlyuVzoOA0qKSnB6NGjsWzZMri7uwsdRytqtRoeHh5YunQpevXqhfDwcMycObPRo7AKKTk5GXPmzMHixYtx+vRp/P3339ixYwdmzZoldDSjEPzKyVjDsPBNl9w15s2bh9jYWOzduxddu3Y1ZMw6aZv9ypUruHbtGoYPH65ZplarAQASiQTZ2dnw9/c3bGjo9jP38vKCjY0NxGKxZtkzzzyDvLw8VFVVQSqVGjRzDV2yf/XVVxg9ejTee+89AECXLl1QWlqKcePGYebMmbwPUcKX+s5PJyenRl81ASZw5WSuw7DokhsA5s6di1mzZiExMREBAQHGiPoEbbN36NABGRkZSE9P17QRI0Zg4MCBSE9PN9oIpbr8zPv27YvLly9riikAXLx4EV5eXkYrTIBu2cvKyp4oQDVFlpnwI7G8nZ/a3as3jISEBCaTydiKFStYZmYmGzduHHNxcWF5eXmMMcZGjx7Npk2bpln/yJEjTCKRsHnz5rGsrCwWHR0tWFcCbXLHxsYyqVTKNmzYwO7cuaNpJSUlRs2tS/Z/E+rTOm1z5+bmMkdHR/bxxx+z7Oxstn37dubh4cG+//57k88eHR3NHB0d2Zo1a1hOTg7bvXs38/f3Z6NGjTJq7pKSEpaWlsbS0tIYAPbTTz+xtLQ0dv36dcYYY9OmTWOjR4/WrF/TlWDy5MksKyuLxcfHm29XAsYYW7hwIfP19WVSqZQFBgaylJQUzddCQkJYZGRkrfXXrVvH2rVrx6RSKevUqRPbsWOHkRM/ok3uli1bMgBPtOjoaOMHZ9r/zB8nVHFiTPvcR48eZUFBQUwmk7HWrVuz2bNnM6VSaeTUj2iTvbq6mn3zzTfM39+fyeVy5uPjw8aPH88KCwuNmnn//v11vm5rskZGRrKQkJAntunevTuTSqWsdevW7I8//tD6uDRkCiHEJAl+z4kQQupCxYkQYpKoOBFCTBIVJ0KISaLiRAgxSVScCCEmiYoTIcQkUXEihJgkKk6EEJNExYkQYpKoOBFCTBIVJ0KISfp/c6lj5IMliDMAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 300x300 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAScAAAEYCAYAAAAedjA5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA3v0lEQVR4nO3deVxU9d4H8M+cMwsgqyCbggiKmmsuIK6PStrVq3m75ZIp2WIm3rpi5tJCqYmPadk1vJWltqhYlpbL44YauZBXhATBFZRFQUkZEHSYOfN7/vBhnlBAZubMnFm+79fr91IOZ/kwzPnyO2fO+R0ZY4yBEEJsDCd1AEIIaQgVJ0KITaLiRAixSVScCCE2iYoTIcQmUXEihNgkKk6EEJtExYkQYpOoOBFCbBIVJ0KITTK6OKWlpWHMmDEIDg6GTCbD9u3bH7rM4cOH0atXL6hUKrRv3x4bNmwwISohxJkYXZyqq6vRo0cPJCcnN2v+goICjB49GkOHDkVWVhb++c9/4sUXX8TevXuNDksIcR4yc278lclk2LZtG8aNG9foPPPmzcOuXbuQk5NjmDZx4kRUVFRgz549pm6aEOLg5JbewPHjxxEbG1tv2siRI/HPf/6z0WU0Gg00Go3ha71ej5s3b8LX1xcymcxSUQkhJmCMoaqqCsHBweA48U5jW7w4lZaWIiAgoN60gIAAVFZW4s6dO3B1dX1gmaSkJLz33nuWjkYIEVFRURHatGkj2vosXpxMsWDBAiQkJBi+VqvVCA0NRVFRETw9PSVMRgi5X2VlJUJCQuDh4SHqei1enAIDA1FWVlZvWllZGTw9PRvsNQGASqWCSqV6YLqnpycVJ0JslNinXCx+nVNMTAxSU1PrTdu/fz9iYmIsvWlCiB0zujjdvn0bWVlZyMrKAnDvUoGsrCwUFhYCuHdINnXqVMP8M2bMQH5+Pt544w2cPXsWa9aswXfffYfZs2eL8xMQQhwTM9KhQ4cYgAdaXFwcY4yxuLg4NmTIkAeW6dmzJ1MqlSw8PJytX7/eqG2q1WoGgKnVamPjEkIszFL7p1nXOVlLZWUlvLy8oFar6ZwTITbGUvsn3VtHCLFJVJwIITaJihMhxCbZ5EWYxH7k5eVh+/btOHr0KM6cOYPCwkLo9XoA96574Xm+3vUv95/ivP97giAY5uF5HqGhoejatSsGDx6McePGoX379lb4qYgtoBPi5KH0ej127dqFzz77DAcPHsSdO3cMRUen0xnm4zgO3t7e8Pf3h7+/P/z8/ODl5QV3d3e4ublBpVJBoVCA4zgwxqDX66HVaqHRaFBTU4Pbt29DrVajvLwc169fx/Xr11FRUVGvoMnlckMRc3d3R2xsLGbOnInY2Fi671Iilto/qedEHnDmzBm8/fbb2LZtG2Qy2QO9HQAQBKHe1+7u7mjVqhWUSiVUKhV0Oh1qamrg4uICd3d3eHh4ICgoCKGhoQgLC0Pbtm2hUCgemkWr1aK0tBSFhYW4fPkyLl26hAsXLiAnJwdnz57F9u3bsX37dsjlcgiCAJlMhvHjx2PJkiWIiIgQ7TUh1kc9J4LLly9j2rRpOHz48EPnVSqV0Ov19XpMTeF5HjzPP7AMx3EICQlB9+7d0aNHD/Ts2RN9+vRBaGhos3tAer0ely9fRlZWFk6dOoUTJ04gPT0dVVVVkMlk4DgOgiDg6aefxqpVqxAcHNys9RLjWGr/pOLkpJKTkzFr1iypY4DjOPA8D61WCwDw9fXFwIEDMXjwYAwdOhQ9evQwahgOvV6Ps2fP4tdff8Uvv/yC/fv3o7y83LCOwMBAfPrppxgzZoxFfh5nZLH9U9RLOi2ErhAXx+uvv97g1f221nieZxzHMQDMy8uLTZgwgX399desvLzc6J9Zr9ez7Oxs9sEHH7DBgwcznucN22jRogVLSUmxwCvtXCy1f1JxcnCbNm2SvNiY2xQKBQPAOI5jQ4YMYWvWrGFlZWUmvR63bt1iGzduZOPGjWNKpdJQqNq0acOysrJEfvWdAxUnKk7NdvfuXebv7y95UbFUr0omkzGO49iIESPYpk2bWE1NjUmvU2VlJfv222/ZiBEjDOsEwKZNm8Z0Op3IvxXHRcWJitND7d27V/LiYc0ml8sZAObp6cleffVVlpuba/Jrd/XqVbZ8+XIWHh5uWHdgYCC7ePGiiL8hx0TFiYpToxYuXCh5oZC61RWq//qv/2I///wzEwTBpNdSr9ezgwcPsr/97W+M4zjDOaqdO3eK/FtzHFScqDg94JlnnpG8KNhaqytSERERbO3atezu3bsmv75Xrlxhc+fOZS1atGAymYwBYOvWrRPxN+gYqDhRcTKYPHmy5EXA1ltdMQkMDGTJyclmFSm1Ws2WL1/OfH19qUg1gIoTFSf2zjvvSL7T21uTyWRMJpOxoKAg9sUXXzCtVmvy619TU8M+/vhj5ufnZyhSe/bsEfE3bJ+oODlxcXK2E92WKlIAWGRkJNuxYwfT6/Um/z5qamrYBx98wLy8vAyf8F26dEnE37h9oeLkhMWpqqpK8p3a0VrdCe5hw4axnJwcs34/t27dYvPmzWMKhYLJ5XLWvXt3p7wEwVL7J43nZKMGDRok+nPAyP/fsJyWlobu3btjzpw5qKqqMmld3t7eWLZsGc6fP4+xY8fi9OnTUKlUWL58uZiRnZeopc5CnKnndOTIEcl7F87UOI5jgYGBbPv27Wb/7g4cOMDat29vOIQsLi4W4R1h+6jn5AT8/f0xcOBAqWM4Fb1ej+vXr2PcuHF4+umncePGDZPXNXz4cOTk5OC9996DQqFAWFgYpk+fLmJaJyNqqbMQR+85HTt2TPIeBLV756N8fHzYDz/8YPbvNC8vj/Xr18+w7pKSEhHeKbaJek4OqlevXujfv7/UMQjunY9Sq9X4+9//jueee87kc1EA0KlTJxw5cgQrV66EQqFA27ZtsWjRIhHTOj4az0kiNTU1aNGihdQxSCN4nkdISAi2bt2K3r17m7WuM2fOYMKECcjLy4O3tzeuX78OnudFSio9em6dA1m/fj0VJhsnCAKKi4vRr18/rF69usGhipurS5cuyMjIwKuvvoqbN29CoVAgKytLvLCOStSDRAtxpHNOPXr0kPzcCjXj24QJE1hVVZXZv/+ffvqJeXp6MrlczhITE81/Q9kAehy5AxzW0dNB7BfHcejYsSN27Nhh9oMTCgoKMG7cOOTk5KBHjx44deqUSCmlQYd1dqy8vJwKk53T6/W4cOECevfujUOHDpm1rnbt2iE9PR3PPPMMMjMzIZPJUFtbK1JSx0HFycIOHDiAVq1aSR2DiECn0+H27dt47LHH8OWXX5q1LldXV3z99ddYsWIFZDIZXFxcUFhYKFJSx0DFyYKWLl2Kxx57TOoYRESCIEAQBLz44ot46623zDpRLpPJMGfOHPz8889QqVQIDw9HWlqaiGntnKhnsCzEHk+IT5w4UfKTuNQs36ZOncpqa2vNfr+cPHmS+fr6Mrlczr766isR3oHWQ6MS2FFxio6OlnynoWadJpPJ2KhRo0x+yMKfXbp0ibVr147xPM+SkpJEeCdaBxUnOylOdQPkU3OexnEcGzx4MKusrDT7/VNWVsa6devGeJ5n8+bNE+EdaXl0+4odCA4ORn5+vtQxiJXp9XocPXoUw4cPh1qtNmtd/v7+SEtLQ+/evbFixQq89tprIqW0P1ScRNKmTRtcu3ZN6hhEIoIg4NSpU6IUKG9vbxw4cADR0dFITk522gJlUnFKTk5GWFgYXFxcEB0djRMnTjQ5/6pVq9CxY0e4uroiJCQEs2fPxt27d00KbIsiIyNRUlIidQwiMUEQkJWVhdjYWFRWVpq1Lg8PD+zdu9dQoBYuXChSSjti7HFgSkoKUyqVbN26dezMmTPspZdeYt7e3o0+Hnrjxo1MpVKxjRs3soKCArZ3714WFBTEZs+e3ext2vI5p4EDB0p+zoOabTWe59mAAQNYdXW12e+vyspK1qdPH8bzPPvwww9FeMeKz2ZOiEdFRbH4+HjD14IgsODg4EY/XYiPj2fDhg2rNy0hIYENGDCg2du01eIUFxcn+Y5AzTYbz/Ns5MiRTKPRmP0+u3nzJuvcuTPjeZ5t2bJFhHeuuGzihHhtbS0yMjIQGxtrmMZxHGJjY3H8+PEGl+nfvz8yMjIMh375+fnYvXs3Ro0a1eh2NBoNKisr6zVb869//QtfffWV1DGIjRIEAfv378dzzz0HvV5v1rp8fHxw8OBBBAcHY/LkyQ89jeIwjKlkJSUlDAA7duxYvelz585lUVFRjS738ccfG55QAYDNmDGjye0kJiY2+NfIVnpONHIlNWPanDlzRHnfnT9/nvn4+DCe59m1a9dEWacYbKLnZIrDhw9j6dKlWLNmDU6dOoUff/wRu3btwuLFixtdZsGCBVCr1YZWVFRk6ZjNVllZSSNXEqOsXLkSq1evNns9HTp0wO7du8HzPIKDgw1PknFYxlQyjUbDeJ5n27Ztqzd96tSpbOzYsQ0uM3DgQPb666/Xm/bNN98wV1dXJghCs7ZrS+ecYAN/ianZX5PJZGznzp2ivAe3bNnCADC5XC7K+sxlEz0npVKJ3r17IzU11TBNr9cjNTUVMTExDS5TU1MDjqu/mbohSpntDyVVDw17Qkwlk8nw9NNPIycnx+x1jR8/Hu+99x50Oh2GDx8uQjobZWw1S0lJYSqVim3YsIHl5uay6dOnM29vb1ZaWsoYY2zKlCls/vz5hvkTExOZh4cH27x5M8vPz2f79u1jERERbPz48c3epi30nIYPHy75X19q9t3kcjkLDQ1l5eXlZr8fBUFg48aNYzzPs+TkZBHe4aazmUsJGGNs9erVLDQ0lCmVShYVFcXS09MN3xsyZAiLi4szfK3Vatm7777LIiIimIuLCwsJCWEzZ85kt27davb2pC5OO3bskPyNTc0xGs/zbPjw4aI8tryyspK1b9+eyeVylp+fL8I73TQ0TK9Ew/TqdDooFAqrbpM4NplMhvnz52Pp0qVmrysvLw+9evWCRqMx+5IFU9EwvRKhwkTExhhDUlISdu7cafa6OnfujLVr14Ixhq5du4qQznZQcWrCX//6V6kjEAfFcRyeffZZUS6TefbZZzF16lTk5eVh/fr1IqSzDXRY14j8/Hyzn7JBSFPkcjn69u2LtLQ0yOVys9ZVVVWFbt26obi4GGq12qrPRaTDOiujwkQsTafTIT09He+//77Z6/Lw8MDmzZuh1+vRsmVLEdJJj4pTA8x9/DQhzcUYw6JFi/Dbb7+Zva6YmBjMmzcPWq0WK1asECGdtOiw7j6XL19Gu3btLLoNQv5MLpejbdu2OH36NNzc3Mxal0ajQbdu3ZCfn4+qqiq4urqKlLJxdFhnJVSYiLXpdDoUFBTgzTffNHtdKpUKGzZsgF6vR5s2bURIJx0qTn/y8ssvSx2BOCm9Xo+PP/4YR48eNXtd/fv3xyuvvAK1Wo0dO3aIkE4adFj3f/R6veGeP0KkwPM8wsPDkZ2dDZVKZda61Go1IiIicPPmTYtfnEmHdRbm4uIidQTi5ARBwMWLF0W5ctzLywurVq0CYwxTp04VIZ31Uc8JQEFBAcLDw0VfLyGmkMvlOHPmDCIjI81aD2MMMTExOHnypEVPjlPPyYKoMBFbM3PmTLOHFJLJZFi9ejUEQUCXLl1ESmY9Tl+cfvrpJ6kjEFKPTqdDamoqtm3bZva6+vbti0mTJqGoqAhXrlwRIZ31OP1hHQ0gR2wRx3EIDg7GhQsXzD4fWlBQgMjISHAcB41GI1LC/0eHdRawZs0aqSMQ0iC9Xo+SkhKsWrXK7HW1a9cOM2bMgCAIOH36tPnhrMSpe07UayK2zs3NDQUFBfD39zdrPdeuXUNYWBhqa2tFHx6bek4i+/zzz6WOQMhDaTSaJp9U1FxBQUGYOXMmeJ4XZRxza3DanhP1moi9kMvlOHfunNmfKl+9ehVhYWGQyWSinnuinpOI9u/fL3UEQozy7rvvmr2O4OBgPP/889Dr9bh69ar5oSzMKXtO1Gsi9kYmkyEvLw8dO3Y0az0XL15EZGQk2rZti4KCAlGyUc9JJOXl5VJHIMRoPM+Lcu6pffv2GDt2LIqLi23+icFOV5zM/dSDECnodDps2rQJFy9eNHtds2fPhk6nw6RJk0RIZjlOV5zs4CiWkAbxPI/ly5ebvZ7BgwfjkUcewdatW0VIZTlOVZzs9e5sQoB7vaf169ebfTJbJpPhH//4BwDY9HhPTlWcvvnmG6kjEGIWxhhWr15t9nqeeeYZqFQqjBs3zvxQFuI0xam0tFTqCISYTRAEJCcno7q62qz1eHp6YuLEieA4DlqtVqR04nKa4kTDohBHcfv2bVGOAqZNmwadTodXXnlFhFTic5rrnOjaJuIoZDIZIiMjkZeXZ9b7Wq/Xo23btiguLjbrgyK6zskMYgwaT4itYIzh3Llz+OWXX8xaD8dxmDJlCniex61bt0RKJx6nKE4DBw6UOgIhopLL5fj000/NXs/EiRMhCIJNfpLtFId1dEhHHJFcLse1a9fg5+dn8joYY2jfvj3y8/NNPrSjwzoT5eXlSR2BEIsQBAHffvutWeuQyWQYP3485HI57ty5I1IycTh8cYqOjpY6AiEWs27dOrPX8cQTT0Cn0+Gf//yn+YFE5PCHdXRIRxzd77//ju7du5u8vF6vh5+fH27dumXSoR0d1pnA0k86JURqcrkcGzduNGsdHMfhr3/9K+RyuUipxGFScUpOTkZYWBhcXFwQHR2NEydONDl/RUUF4uPjERQUBJVKhcjISOzevdukwMaIj4+3+DYIkZJOp8O3335r9g3tI0aMgE6nw6FDh0RKJgJmpJSUFKZUKtm6devYmTNn2EsvvcS8vb1ZWVlZg/NrNBrWp08fNmrUKHbkyBFWUFDADh8+zLKyspq9TbVazQAwtVptVFYA1Kg5RTt+/LhR+8b9rl69ygCwXr16Gb2sqfvnwxhdnKKiolh8fLzha0EQWHBwMEtKSmpw/n//+98sPDyc1dbWmhySihM1ao03uVzO5syZY/L+VSciIoIBRpcEixUnow7ramtrkZGRgdjYWMM0juMQGxuL48ePN7jMzz//jJiYGMTHxyMgIABdu3bF0qVLLT4KH51vIs5Cp9Phu+++M/vQLjY2FgqFQqRU5jOqOJWXl0MQBAQEBNSbHhAQ0Ohd//n5+di6dSsEQcDu3bvx9ttvY+XKlViyZEmj29FoNKisrKzXjLVw4UKjlyHEXhUVFSE7O9usdQwcOBBarRa//fabSKnMZEw3q6SkhAFgx44dqzd97ty5LCoqqsFlOnTowEJCQphOpzNMW7lyJQsMDGx0O4mJiQ12X43pNja0PDVqjtp4nmdLlixp9v7RkIsXLzIAbPTo0UYtZxOHdX5+fuB5HmVlZfWml5WVITAwsMFlgoKCEBkZCZ7nDdM6d+6M0tJS1NbWNrjMggULoFarDa2oqMiYmIQ4Hb1ej59//tmsdYSHh8PLywu7du0SKZV5jCpOSqUSvXv3RmpqqmGaXq9HamoqYmJiGlxmwIABuHjxYr1zQOfPn0dQUBCUSmWDy6hUKnh6etZrhJDGMcbwn//8x6zRBWQyGaKiomznwmVju1opKSlMpVKxDRs2sNzcXDZ9+nTm7e3NSktLGWOMTZkyhc2fP98wf2FhIfPw8GCzZs1i586dYzt37mT+/v5GdUGN7TbevXtX8m42NWpStO+//964Hfo+8+fPZ3K53KhlbOZSAsYYW716NQsNDWVKpZJFRUWx9PR0w/eGDBnC4uLi6s1/7NgxFh0dzVQqFQsPD2fvv/9+vXNQD2PsDz9nzhzJ3yTUqFm7KRQK9vLLLzd7v2rI5s2bGYB6+/TD2FRxsjZjf3hXV1fJ3yjUqEnR2rVrZ9a+lpOTwwCw559/vtnL2MQJcXtha0M/EGItBQUFZj06qkOHDuA4Dvv37xcxlWkcsjgR4szS0tJMXlapVCI0NNQmPiGn4kSIA1EoFDh27JhZ6+jUqZNIacxDxYkQB6LVanHkyBGz1hEREWETt7E4XHG6fv261BEIkVR2djbu3r1r8vJhYWHQ6/WS35/qcMUpOTlZ6giESEqn0+H06dMmLx8SEgJBEJCbmytiKuM5XHH66aefpI5AiKRkMhlOnTpl8vKtW7cGgHp3gkjB4YqT1NWeEKnJ5XJkZmaavHzdfbIZGRliRTKJwxUnrVYrdQRCJKXVas0qTq1atQIAnDt3TqxIJnG44kQIAc6cOWPy4HOenp7gOE7ya52oOBHigGpqalBSUmLSsjKZDJ6enigvLxc5lXGoOBHioEw9LKusrERFRYXkp0ioOBHigGQyGc6fP9/s+XU6HU6ePIkPP/wQHTt2tGCy5rOtp+gRQkShUChw8eLFBr+n1+tx5coV5ObmIisrC7/88guOHDmCO3fugONsp79CxYkQB6TT6ZCVlYXMzEzk5+fj7NmzyMvLw+nTp3H+/HloNBoAAM/zYIwZrgaX+qrwP6PiRIiDUSgU0Gq1OHjwIHr16gUAhkeN63S6evNa+hFt5qDiRIgd4XkePM9Dp9PV6+XI5XK0bt0a7du3R0REBL799lvU1NQYvn9/UbIHVJwIsQFyuRwcx0EQhAd6M0qlEv7+/ggNDUXbtm3Rpk0btGnTBiEhIYb/BwQE1DtfFBUVhRdffNHaP4aoqDgRYgEcx0Eul4MxBp1O98AFkRzHwdvbGwEBAQgKCkJQUBACAgIQGBiIoKAgBAcHG6Z7eXkZ/USUoKAgMX8cSVBxIuQh6no1jDEIgtDgSWO5XA5vb2/4+fnB398fAQEBaNWqFXx9feHv7w9/f3+0atXK8K+vr2+9ZzmK7f6nctsjKk7E4clkMsjlcshksiYLDHCvyHh6esLb2xu+vr7w8/ODn58fWrZsaWi+vr6G/9d9393d3Xae94Z7D8C1d1SciM3iOA48zxuKSt0AaE3dM6ZQKODh4WEoMC1btoSPjw+8vLzg7e0Nb29veHl5Gab5+PjA29sbPj4+8PHxgZubm00VGVO1bNlS6ghmc9ji5ObmhpqaGsOnG409+pyYj+d5cBxn2KnrCklzP6Z2cXGBm5sbWrRoYSgsXl5e8PT0hIeHh6H9+ev7nwhd11QqlSV/VLvh7u4OjuNs6rolYzlkcQoICMClS5ewf/9+XLhwAQUFBbh06RIuXryIoqKiRu8Z4jjOsJP9eUerOxQw9S5vsdXl+3NBuP+vfV3uukJhanaFQgGVSgUXFxe4uLjA1dXVUEjc3Nzg7u4ONzc3uLq6Gqb9+d+Gmru7e71/bemqZEchk8nQokULVFVVSR3FZA5XnLy9vaFWq9GiRQuMGzfuge8zxlBdXY2KigpDu3XrluH/NTU1uHPnDu7evYs7d+7U+//du3eh0+mg1WrrNZ1Oh9ra2nqF4P5/7y8idf/W9TrqWl1Pj+d5yOVyw79yuRwKhQIKheKB/yuVSsPXdcWk7l+lUgmlUmn4v0qlqtfqio6Li8sDRcjFxYUKhx1zc3Oj4mRL2rRpg5ycHNy5cweurq4PfF8mk8Hd3R3u7u5o06aNBAkJsY6G3v/WXN5cDvdnMTIyEgBQVlYmcRJCpOXi4mLW8o888ohISUzjcMWpT58+AIBr165JnIQQaSmVSrOWHz16tEhJTONwxWnEiBEAgOLiYomTECItcy/ynDFjhkhJTONwxenRRx8Fx3EoLCyUOgohdk3qW2AcrjjVfeJVUFAgdRRCJGXP1zgBDlicgHuPxrlw4YLUMQiRlC2P1dQcDlmcACAvL0/qCIRIqm60S3vlkMWpXbt2KCkpwd27d6WOQohk7P3975DFadSoUdDr9ZI/sZQQKf15JExjPfXUUyImMY1DFqf4+HgAwOnTpyVOQoh0qqurTV72X//6l4hJTGNScUpOTkZYWBhcXFwQHR2NEydONGu5lJQUyGSyBu95E1Pnzp2hUCiQlZVl0e0QYqu0Wq1Zh3VSX0YAmFCctmzZgoSEBCQmJuLUqVPo0aMHRo4cievXrze53OXLl/H6669j0KBBJoc1hlarxX/+8x+rbIsQW1NRUSF1BLMZXZw+/PBDvPTSS5g2bRoeeeQRfPrpp3Bzc8O6desaXUYQBEyePBnvvfcewsPDzQrcXC4uLsjIyLD7j1MJMUV5ebnUEcxmVHGqra1FRkYGYmNj/38FHIfY2FgcP3680eUWLVoEf39/vPDCC83ajkajQWVlZb1mrAkTJqCmpga5ublGL0uIvXvYkUxThg0bJmIS0xlVnMrLyyEIwgODpwcEBKC0tLTBZY4cOYIvv/wSa9eubfZ2kpKS4OXlZWghISHGxAQAJCYmQiaT4dixY0YvS4i9a2x/bI7vv/9exCSms+indVVVVZgyZQrWrl1r1IDrCxYsgFqtNrSioiKjt92uXTtwHIe0tDSjlyXE3pWUlJg8UKCtjD9u1GBzfn5+4Hn+gbGSysrKEBgY+MD8ly5dwuXLlzFmzBjDtLr7feRyOc6dO4eIiIgHlqsbpdFcgiDgwIED9UaiJMQZFBUVQS6X2/XY+UaVVqVSid69eyM1NdUwTa/XIzU1FTExMQ/M36lTJ2RnZyMrK8vQxo4di6FDhyIrK8ukwzVjjBkzBtevX8f58+ctuh1CbM3ly5cbHSu/KV26dLFAGtMYPUxvQkIC4uLi0KdPH0RFRWHVqlWorq7GtGnTAABTp05F69atkZSUBBcXF3Tt2rXe8t7e3gDwwHRL+Oijj7Bz507s378fHTt2tPj2CLEVFy5cMOmhFvv377dAGtMYXZwmTJiAGzdu4J133kFpaSl69uyJPXv2GE6SFxYW2syg+BEREZDJZNi9ezdmzZoldRxCrIIxhkuXLpm0rC1cfFlHxmzleUdNqKyshJeXF9RqNTw9PY1a1t3dHTqdDrdu3ZJ8wHZCrKG4uNjkUyamlANz9s+m2EYXx4KWLl0KjUZT7zwZIY7M1OGCPvjgA5GTmMfhi9OsWbMgl8uxfft2qaMQYhU5OTkmnVp5/fXXLZDGdA5fnDiOg06nw9atW0369IIQe3P69GmzH25gCxy+OAHAm2++CbVajcOHD0sdhRCLO3nypNF/iGfOnGmhNKZz+BPiwL1rsVQqFZ599lmsX7/eAgkJsQ13796Fu7u70Te8m1MG6IS4GeoO7b7//nvcuXNH6jiEWExmZqbDjMThFMUJANatW4fq6mo6MU4cWnp6utEnw5OTky2UxjxOU5ymTZsGnueNGh2BEHvz66+/Gn0fqS2ebwKcqDgBQPv27XHo0CHk5+dLHYUQ0THGcPjwYaMO6xQKhQUTmcepitO+ffuo90Qc1pkzZ3Dr1i2jlikuLrZQGvM5VXEKDQ2FIAj49NNP7f6ZXoTcLzU11ehDOn9/fwulMZ9TFScA+Pzzz1FRUYHNmzdLHYUQUe3Zs8eo4rRixQoLpjGfU1zndD+O49C5c2fk5OTQIHTEIdy5cwc+Pj5GPYJcrF2frnMS0eTJk5Gbm2tTY9cQYo6DBw8aVZgGDx5swTTicMri9NVXX4Hnebz//vtSRyFEFD/99JNRn7z98ssvFkwjDqcsThzHYcCAAUhLS6OnsxC7JwgCfvjhh2bfT2dLA8o1xSmLE3BvOFK5XI7ExESpoxBill9//RU3b95s9vy2fPnAnzltcap7WMOBAweo90Ts2ubNmyGXN2/EbV9fX5sZRvthnPLTujq1tbVwc3NDTEwM0tLS6JM7Ync0Gg38/f2b/VRsQRBEL070aZ0FKJVKjBgxAkeOHMGePXukjkOI0Xbu3NnswhQZGWk3vSbAyYsTcO+Xy3EcEhISHGaoCeI8vvjii2Yf0p07d87CacTl9MWJ4zi8/vrrOHv2LL788kup4xDSbEVFRdi7dy90Ot1D5500aZIVEonLqc85/ZlMJoOPjw/y8/MND/4kxJa98847WLp0abN6/Jbczemck4Xt3r0blZWVeOedd6SOQshD1dbWYs2aNc0qTD/88IMVEomPitP/+ctf/gI3Nzd88sknOHXqlNRxCGnSli1b8McffzRr3ieffNLCaSyDitOfFBcXg+M4vPDCC3RynNgsxhiWL1/erE/eampqrJDIMqg4/YmnpydmzZqFrKwsfPzxx1LHIaRB+/btQ05ODvR6fZPzPfHEE3B1dbVSKvHRCfEGyGQyqFQq5OTkoH379hbfHiHGGDhwINLT0x/au7fWrk0nxK0oPz8fgiBgypQpdHhHbMovv/yCo0ePPvR9WVBQYKVElkPFqQHt2rVDXFwc0tPTsXLlSqnjEALgXk9o4cKFD73ocujQoQgLC7NOKAuiw7omyGQyyOVy/Pbbb+jVq5fVtktIQ/bs2YO//OUvD53P2rs0HdZJoLy8HIwxPP3007h9+7bUcYgTEwQBc+fOBc/zTc5XUVFhnUBWQMWpCb6+vkhOTsbly5fx8ssvW/0vEiF1vv76a+Tk5DR5rmnRokXw8vKyYirLosO6Znj00UeRlZWFzz77DNOnT7f69olzq6qqQnh4OP74449G/0C6ubmhurraysnuocM6CWVmZgIA4uPjceLECYnTEGezaNEi3Lx5s8meu1SFyZJMKk7JyckICwuDi4sLoqOjm9xh165di0GDBsHHxwc+Pj6IjY21yx38xo0bYIxh7NixKC0tlToOcRLZ2dn46KOPmrzg0pgheu2J0cVpy5YtSEhIQGJiIk6dOoUePXpg5MiRuH79eoPzHz58GJMmTcKhQ4dw/PhxhISEYMSIESgpKTE7vDX5+flhx44dKC8vx9ixY+mJwcTiBEHACy+80OQIrSkpKfDx8bFiKitiRoqKimLx8fGGrwVBYMHBwSwpKalZy+t0Oubh4cG++uqrZm9TrVYzAEytVhsbV3Rz5sxhHMexiRMnMr1eL3Uc4sA++ugjBqDR9uSTT0odkTFmuf3TqJ5TbW0tMjIyEBsba5jGcRxiY2Nx/PjxZq2jpqYGWq0WLVu2bHQejUaDysrKes1WrFixAn379kVKSgoNr0Is5vz585g/f36j3/fy8rLboVCay6jiVF5eDkEQEBAQUG96QEBAs8/DzJs3D8HBwfUK3P2SkpLg5eVlaCEhIcbEtLj09HS4urpiyZIl+Pzzz6WOQxyMVqvF5MmTm7xswJGuZ2qMVT+tW7ZsGVJSUrBt2za4uLg0Ot+CBQugVqsNraioyIopm6fuoswZM2bgxx9/lDgNcSSLFy9GRkZGo8PvMtu/+kcURhUnPz8/8DyPsrKyetPLysoQGBjY5LIrVqzAsmXLsG/fPnTv3r3JeVUqFTw9Pes1W8NxnOHZ9BMmTMC+ffskTkQcwaFDh7BkyZJGC1Bzn+rrCIwqTnUPokxNTTVM0+v1SE1NRUxMTKPLLV++HIsXL8aePXvQp08f09PaGKVSiT/++AOCIGDs2LF28fx5YruuXbuG8ePHNzqIXGVlZbOftOIQjD2DnpKSwlQqFduwYQPLzc1l06dPZ97e3qy0tJQxxtiUKVPY/PnzDfMvW7aMKZVKtnXrVnbt2jVDq6qqavY2benTuoaUlJQwjuOYi4sLS0tLkzoOsUMajYbFxMQwuVze4CdzV69elTpioyy1fxpdnBhjbPXq1Sw0NJQplUoWFRXF0tPTDd8bMmQIi4uLM3zdtm3bBl/sxMTEZm/P1osTY4xduXLFUKAOHTokdRxiZ15++WUmk8ka3FeuXLkidbwmWWr/pHvrRFRUVISwsDDI5XL89NNPePzxx6WOROzA6tWr8eqrrzb4vStXriA0NNTKiYxD99bZgZCQEFy9ehVarRZjxozBli1bpI5EbNyOHTvw2muvNfi90tJSmy9MlkTFSWQBAQGoqKiATqfDpEmTsHr1aqkjERt1/PhxPP300w1+r6qq6oHrCZ0NFScL8PT0RG1tLRhjePXVVzF37tyHPimDOJfs7Gw8/vjj0Ol0D1w2IAgC3N3dJUpmO6g4WYhCoQBjDL6+vlixYgWefPJJhxzWghjv3LlzGDp0KGpqah64Cpwx1qzn0TkDehUsrLy8HLGxsdixYwf69euHwsJCqSMRCZ07dw6DBg2CWq2udwV4UFCQ01z53VxUnKxg//79WLx4MfLy8tCzZ08cPnxY6khEAjk5ORgwYABu3bpVrzBNmzYNV69elTCZbaLiZCULFy7Er7/+CrVajeHDh+ODDz6gv5RO5LfffsPAgQMNH5bU2bt3L9atWydhMttF1zlZWXV1teFk56hRo/D111/D19dX4lTEknbt2oWnnnrqgQEKa2pq7Ppx4XXoOicH0aJFCzDGMGrUKOzduxddu3alwzwHlpycjDFjxtQrTG3atAFjzCEKkyVRcZLIrl27sG3bNly/fh3Dhg3DG2+8YRjlgNg/rVaLmTNnYtasWfUO37/44gubHALIFtFhncTu3r2Lli1b4u7du+jYsSO++eYbhxq5wRldu3YNf//73x8YHba2thYKhUKiVJZDh3UOysXFBTU1NViyZAkuXLiA6OhozJs3DzU1NVJHIyZITU1FcHBwvcI0bdo0MMYcsjBZEvWcbIharUarVq0gCAJCQkLw2WefYeTIkVLHIs1QW1uLxMRELFu2rN70yspKeHh4SJTKOqjn5AS8vLxQW1uLTz75BEVFRXj88cfx5JNP4sqVK1JHI03Izs6GSqWqV5jmzZsHxpjDFyZLop6TjRIEAX369EF2djZ4nse8efPwxhtv0D1XNkSj0WDGjBnYsGGDYZpcLodGo3GqW1Co5+RkeJ5HZmYm8vLyUFtbi/fffx/t2rXDZ5991ujA98R66h7S8efClJ2dDa1W61SFyZLoVbRxHTp0AGMMmzZtQnl5OWbMmIEOHTpg06ZNTT46iFhGTk4OZDIZnnzyScO0FStWgDGGrl27SpjM8VBxshMTJkwAYwyLFi3ClStXMHnyZHTq1AkbN26knpQV5OfnQyaToVu3boZpEydOBGMMc+bMkTCZ46LiZGfefvtt6PV6zJkzB/n5+Xj22WcRERGBNWvW0OUHFpCZmQmZTIaIiAjDtMGDB4Mxhs2bN0uYzPFRcbJTK1asgCAIePfdd1FUVIRZs2YhODgYCxcuRHFxsdTx7N5nn30GmUyGXr16GaYNGzYMjDF6BJi1iPq4BAuxh6evSO3rr79mABjP84zjOPbEE0+wPXv2MEEQpI5mN6qrq1lwcPADTz957bXXpI5m0yy1f1LPyUFMmTIFjDHk5uYiICAAu3btwuOPP46QkBC88847uHjxotQRbVZeXh4WLFiAFi1a1BtXKS0tDYwxrFq1SrpwToyuc3JQer0eCQkJ+Pjjj8HzPARBQO/evfHss8/iqaeeQps2baSOKKmCggJs3boV33zzjeFaMkEQ0KpVK1y/fl3qeHbFUvsnFScnUFxcjKeeegonTpyATCaDXq9Hr1698NRTT2HMmDHo0qULZDKZ1DEtijGGrKws7Ny5E99//z2ys7PBcRz0ej0GDx6MrVu3olWrVlLHtEtUnKg4iSInJweTJ0/G6dOnDTtnYGAgRo8ejcceewxDhw6Fv7+/1DFFce3aNRw8eBAHDhzArl27cOPGDUMPqW/fvti8eXO9T+GIaag4UXESXVlZGV555RVs27YNcrnccL1U+/btMWzYMMTExKBfv36IjIy0+aueBUHA2bNnkZ6ejuPHj+PgwYMoKCgAAMPPNmnSJHzyySdo2bKlxGkdCxUnKk4W9z//8z9YsGABfv/993rFys3NDY8++ih69eqFbt264ZFHHkHHjh3h6+tr9cNBxhhu3LiBc+fO4cyZM8jJycHJkyfx+++/G0abrMvet29f/Pd//zeGDh1q1YzOhooTFSerO336NJKSkrBt2zZoNBooFApotVrD9z08PNCuXTt06NABISEhaNOmDYKCguDv7w9fX1/4+PjAy8sL7u7uDx3LqLa2FtXV1aioqMCtW7fwxx9/oKysDNeuXUNJSQkKCwtx4cIFXL58Gbdv3wYAyGQyyOVyaLVatGjRAn/729/w1ltvoWPHjhZ9XUh9VJyoONmE0tJSrF+/Hjt27EBGRgZqa2vBcRx4noder2/0fj+e56FQKCCXyw2HiHq9HjqdDlqttsnlOI4zPBnX1dUVffr0wdixY/H888/TIZoNoOJExcku5Ofn47fffsPvv/+O/Px8lJSU4ObNm6isrERNTQ1qa2shCAI4jgPHcVAqlXB1dYWnpyd8fX3RunVrREREoHv37ujXrx/CwsKk/pHIQ1BxouJEiE2i8ZwIIU6FihMhxCZRcSKE2CQqToQQm2RScUpOTkZYWBhcXFwQHR2NEydONDn/999/j06dOsHFxQXdunXD7t27TQpLCHEeRhenLVu2ICEhAYmJiTh16hR69OiBkSNHNnon97FjxzBp0iS88MILyMzMxLhx4zBu3Djk5OSYHZ4Q4riMvpQgOjoaffv2xSeffALg3oV0ISEh+Mc//oH58+c/MP+ECRNQXV2NnTt3Gqb169cPPXv2xKefftqsbdKlBITYLkvtn3JjZq6trUVGRgYWLFhgmMZxHGJjYx94Lnyd48ePIyEhod60kSNHYvv27Y1uR6PRQKPRGL5Wq9UA7r0IhBDbUrdfin3JpFHFqby8HIIgICAgoN70gIAAnD17tsFlSktLG5y/tLS00e0kJSXhvffee2B6SEiIMXEJIVb0xx9/wMvLS7T1GVWcrGXBggX1elsVFRVo27YtCgsLRf3hLa2yshIhISEoKiqyu8NRe81ur7kB+82uVqsRGhoq+n2ORhUnPz8/8DyPsrKyetPLysoQGBjY4DKBgYFGzQ8AKpUKKpXqgeleXl529Uur4+npaZe5AfvNbq+5AfvNLvaYX0atTalUonfv3khNTTVM0+v1SE1NRUxMTIPLxMTE1JsfAPbv39/o/IQQAphwWJeQkIC4uDj06dMHUVFRWLVqFaqrqzFt2jQAwNSpU9G6dWskJSUBAF577TUMGTIEK1euxOjRo5GSkoKTJ0/i888/F/cnIYQ4FKOL04QJE3Djxg288847KC0tRc+ePbFnzx7DSe/CwsJ63bv+/ftj06ZNeOutt7Bw4UJ06NAB27dvN+q58iqVComJiQ0e6tkye80N2G92e80N2G92S+W2iyFTCCHOh+6tI4TYJCpOhBCbRMWJEGKTqDgRQmySzRQnex2GxZjca9euxaBBg+Dj4wMfHx/ExsY+9Oe0JGNf8zopKSmQyWQYN26cZQM2wtjcFRUViI+PR1BQEFQqFSIjI+3i/QIAq1atQseOHeHq6oqQkBDMnj3b8Hw+a0lLS8OYMWMQHBwMmUzW5H2xdQ4fPoxevXpBpVKhffv22LBhg/EbZjYgJSWFKZVKtm7dOnbmzBn20ksvMW9vb1ZWVtbg/EePHmU8z7Ply5ez3Nxc9tZbbzGFQsGys7NtOvczzzzDkpOTWWZmJsvLy2PPPfcc8/LyYsXFxVbNzZjx2esUFBSw1q1bs0GDBrEnnnjCOmH/xNjcGo2G9enTh40aNYodOXKEFRQUsMOHD7OsrCwrJzc++8aNG5lKpWIbN25kBQUFbO/evSwoKIjNnj3bqrl3797N3nzzTfbjjz8yAGzbtm1Nzp+fn8/c3NxYQkICy83NZatXr2Y8z7M9e/YYtV2bKE5RUVEsPj7e8LUgCCw4OJglJSU1OP/48ePZ6NGj602Ljo5mL7/8skVz3s/Y3PfT6XTMw8ODffXVV5aK2ChTsut0Ota/f3/2xRdfsLi4OEmKk7G5//3vf7Pw8HBWW1trrYiNMjZ7fHw8GzZsWL1pCQkJbMCAARbN2ZTmFKc33niDdenSpd60CRMmsJEjRxq1LckP6+qGYYmNjTVMa84wLH+eH7g3DEtj81uCKbnvV1NTA61Wa/UHQ5qafdGiRfD398cLL7xgjZgPMCX3zz//jJiYGMTHxyMgIABdu3bF0qVLG32Ip6WYkr1///7IyMgwHPrl5+dj9+7dGDVqlFUym0qs/VPyUQmsNQyL2EzJfb958+YhODj4gV+kpZmS/ciRI/jyyy+RlZVlhYQNMyV3fn4+Dh48iMmTJ2P37t24ePEiZs6cCa1Wi8TERGvEBmBa9meeeQbl5eUYOHAgGGPQ6XSYMWMGFi5caI3IJmts/6ysrMSdO3fg6urarPVI3nNyVsuWLUNKSgq2bdsGFxcXqeM0qaqqClOmTMHatWvh5+cndRyj6PV6+Pv74/PPP0fv3r0xYcIEvPnmm80ehVVKhw8fxtKlS7FmzRqcOnUKP/74I3bt2oXFixdLHc0qJO85WWsYFrGZkrvOihUrsGzZMhw4cADdu3e3ZMwGGZv90qVLuHz5MsaMGWOYptfrAQByuRznzp1DRESEZUPDtNc8KCgICoUCPM8bpnXu3BmlpaWora2FUqm0aOY6pmR/++23MWXKFLz44osAgG7duqG6uhrTp0/Hm2++KfoQJWJpbP/09PRsdq8JsIGek70Ow2JKbgBYvnw5Fi9ejD179qBPnz7WiPoAY7N36tQJ2dnZyMrKMrSxY8di6NChyMrKstoIpaa85gMGDMDFixcNxRQAzp8/j6CgIKsVJsC07DU1NQ8UoLoiy2z4lljR9k/jztVbRkpKClOpVGzDhg0sNzeXTZ8+nXl7e7PS0lLGGGNTpkxh8+fPN8x/9OhRJpfL2YoVK1heXh5LTEyU7FICY3IvW7aMKZVKtnXrVnbt2jVDq6qqsmpuU7LfT6pP64zNXVhYyDw8PNisWbPYuXPn2M6dO5m/vz9bsmSJzWdPTExkHh4ebPPmzSw/P5/t27ePRUREsPHjx1s1d1VVFcvMzGSZmZkMAPvwww9ZZmYmu3LlCmOMsfnz57MpU6YY5q+7lGDu3LksLy+PJScn2++lBIwxtnr1ahYaGsqUSiWLiopi6enphu8NGTKExcXF1Zv/u+++Y5GRkUypVLIuXbqwXbt2WTnxPcbkbtu2LQPwQEtMTLR+cGb8a/5nUhUnxozPfezYMRYdHc1UKhULDw9n77//PtPpdFZOfY8x2bVaLXv33XdZREQEc3FxYSEhIWzmzJns1q1bVs186NChBt+3dVnj4uLYkCFDHlimZ8+eTKlUsvDwcLZ+/Xqjt0tDphBCbJLk55wIIaQhVJwIITaJihMhxCZRcSKE2CQqToQQm0TFiRBik6g4EUJsEhUnQohNouJECLFJVJwIITaJihMhxCZRcSKE2KT/BY2aM/en4p5hAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 300x300 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from obspy.imaging.beachball import aux_plane, beachball\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "\n",
    "def plot_beachball(source, size=200, linewidth=1):\n",
    "    \"\"\"\n",
    "    Plot a beachball diagram using given focal mechanism parameters.\n",
    "\n",
    "    :param source: Dictionary containing focal mechanism parameters.\n",
    "    :param size: Size of the beachball diagram.\n",
    "    :param linewidth: Line width of the beachball diagram.\n",
    "    \"\"\"\n",
    "    # Create a figure and axis\n",
    "    fig, ax = plt.subplots(figsize=(3, 3))\n",
    "\n",
    "    # Extract the focal mechanism parameters\n",
    "    strike = source[0]\n",
    "    dip = source[1]\n",
    "    rake = source[2]\n",
    "\n",
    "    # Plot the beachball\n",
    "    beachball([strike, dip, rake], size=size,\n",
    "              linewidth=linewidth, facecolor='black', bgcolor='white', fig=fig)\n",
    "\n",
    "    # Show the plot\n",
    "    plt.show()\n",
    "\n",
    "\n",
    "# 定义原始震源机制\n",
    "strike = 339\n",
    "dip = 15\n",
    "rake = 161\n",
    "\n",
    "\n",
    "# 获取等效节面（辅助面）\n",
    "strike2, dip2, rake2 = aux_plane(strike, dip, rake)\n",
    "\n",
    "\n",
    "# 打印等效解\n",
    "print(\"等效解（辅助面）:\")\n",
    "print(f\"{strike2:.2f}, {dip2:.2f}, {rake2:.2f}\")\n",
    "plot_beachball([strike,dip,rake])\n",
    "plot_beachball([strike2, dip2, rake2])"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "base",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "undefined.undefined.undefined"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
